Extremal black holes in D=4 Gauss-Bonnet gravity 
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We show that four-dimensional Einstein-Maxwell-Dilaton-Gauss-Bonnet gravity admits asymp- 
totically flat black hole solutions with a degenerate event horizon of the Reissner-Nordstrom type 
AdS2 x S 2 . Such black holes exist for the dilaton coupling constant within the interval < a 2 < a%. 
Black holes must be endowed with an electric charge and (possibly) with magnetic charge (dyons) 
but they can not be purely magnetic. Purely electric solutions are constructed numerically and the 
critical dilaton coupling is determined a cr ~ 0.488219703. For each value of the dilaton coupling 
a within this interval and for a fixed value of the Gauss-Bonnet coupling a we have a family of 
black holes parameterized by their electric charge. Relation between the mass, the electric charge 
and the dilaton charge at both ends of the allowed interval of a is reminiscent of the BPS condition 
for dilaton black holes in the Einstein-Maxwell-Dilaton theory. The entropy of the DGB extremal 
black holes is twice the Bekenstein-Hawking entropy. 

PACS numbers: 04.20.Jb, 04.65. +e, 98.80.-k 



I. INTRODUCTION 

String theory suggests higher curvature corrections to the Einstein equations [J 0, Q. Black holes in higher- 
curvature gravity were extensively studied during two past decades [H Q culminating in recent spectacular progress 
in the microscopic string calculations of the black hole entropy (for a review see [g, 0)- ^ n theories with higher 
curvature corrections, classical entropy deviates from the Bekenstein-Hawking value and can be calculated using 
Wald's formalism [H, 0, [ll[ . Remarkably, it still exhibits exact a gree ment with string theory quantum predictions 
at the corresponding level, 

M HI M M SO, [32 



2J, 123J, UM and non-BPS [25|, |2fl 
models with higher curvature terms 

exact classical solutions for static black holes were obtained [20, [U, [24| . Moreover, as was argued by Sen [33, [iy] , 
knowledge of exact global black hole solutions is not necessary to be able to compare classical and quantum results: 
the entropy can be computed locally using the "entropy function" approach based on the typical for supergravities 
attractor property (2f| [U [H, [29|, [3(| • In this case it is tacitly assumed that global asymptotically flat black hole 
solutions exist indeed. Generically, however, the existence of local solutions does not imply possibility to extend them 
to infinity as asymptotically flat black holes. Here we investigate this issue within a simple model of higher curvature 
gravity. 

One of the simplest four-dimensional models with higher-curvature terms is the so-called dilaton-Gauss-Bonnet 
gravity (DGB) which is obtained by adding to the Einstein action the four-dimensional Euler density multiplied 
by the dilaton exponent. As other higher-curvature theories based on topological invariants, this theory does not 
contain higher derivatives and thus is free of ghosts. Black hole solutions in this theory can not be found in analytical 
form, but they were extensively studied perturbatively [U [42j and numerically [43|, [3, EE EE|. More rece ntly 
global properties of DGB black hole solutions were studied using the dynamical system approach [13, [H, IS l50f - 
Stability issues were discussed in [5l|, HH, [H, [H, [55j . In these papers the existence of both neutral and charged 
asymptotically flat solutions with a non-degenerate event horizon and without naked singularities was established. 
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These solutions have the Schwarzschild type event horizon and do not possess an extremal limit. In this respect 
they differs substantially from the dilatonic black holes in the Einstein-Maxwell-dilaton (EMD) theory without the 
GB term [HB, [E3, US, Ull : charged dilatonic black holes do have an extremal limit in which case the event horizon 
shrinks to a point-like singularity. The Bekenstein-Hawking entropy of the extremal dilatonic black holes is zero, while 
quantum theory suggests a non-zero result. The puzzle was solved in several supersymmetric models by showing that 
the horizon of the extremal dilatonic black hole is stretched to a finite radius. In the case of the DGB black holes, 
however, no solution with the degenerate event horizon of finite radius was found so far. 

The aim of the present paper is to study this possibility in more detail. We show that, apart from the known DGB 
black holes with a non-degenerate event horizon, there exist electrically charged solutions with the degenerate horizon 
of the A1IS2 x S 2 type which are asymptotically flat. These new solutions exist only in a limited range of the dilaton 
coupling constant. For other values of this constant, local solutions with the AdS2 x S 2 horizon can not be continued 
to infinity as asymptotically flat black holes: singularity is met in a finite distance outside the horizon. Since the 
DGB theory does not possess S-duality, magnetic solutions differ substantially from the electric ones; in particular, no 
purely magnetic black holes with a degenerate horizon are allowed, though dyonic solutions with a non-zero electric 
charge are possible. 



II. 4D DILATONIC GAUSS-BONNET THEORY 



The action for the four-dimensional Einstein-Maxwell-dilaton theory with an arbitrary dilaton coupling constant a 
modified by the DGB term reads 



S = 7^ / I R - 29 M<W - e 2 <^ (F 2 - a£ GB ) } ^f~g d 4 
where Cqb is the Euler density 

£gb = R 2 — 4i? Atl/ i? A " y + R a fi llv R a ^ llv . 



(1) 



(2) 



This action contains two parameters (we use units G = c = 1): the dilaton coupling a and the GB coupling a. We 
will assume a > 0, a > 0. Solutions for negative a can be obtained changing the sign of the dilaton. Note that the 
Maxwell term is not multiplied by a to facilitate decoupling of the GB term from the EMD action. 

Consider the static spherically symmetric metrics parameterized for further convenience by three functions 
w(r), p(r), a(r): 



ds 2 = -w(r)a 2 (r)dt 2 + 

The scalar curvature and the Euler density then read 

1 



dr 2 
w(r) 



p 2 (r)dQ 2 , 
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The corresponding ansatz for the Maxwell one-form is 

A = —f(r) dt — q m cos 9 dip, 



(3) 

(4) 
(5) 

(6) 



where f(r) is the electrostatic potential and q m is the magnetic charge. Note that the DGB term breaks the discrete 
S-duality which in absence of this term is described by the transformation 



9 M u -» 9^, F - e~ 2a **F, 



(7) 



where F = dA. It is expected therefore that properties of electric or magnetic black holes in this theory will be be 
essentially different. 



A. Reduced action and field equations 



The corresponding one-dimensional Lagrangian is obtained by dropping the total derivative in the dimensionally 
reduced action: 



L = ^[p'(w P y + 1] + a'wpp' - 2aaa- 1 (a 2 wY(wp' 2 - l)0'e 2 ^ - \p 2 ^' 2 + \ P -f 2 ^ ' " 1 
2 2 2 a 



The corresponding equations of motion read: 

8aa [w(wp' 2 - l)0'e 2a *l' - p'(wp)' + 1 - 2wpp" - Aaaw'(wp' 2 - l)0'e 2a0 - wp 2 (j)' 2 
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The Maxwell equation for the form field ([13]) can be directly solved 

/'(r) = g e( rp- 2 e- 2 ^ 

where q e is the electric charge. 



0, 
0, 
0, 
0, 



ft e 2a<p\ = Q_ 



(8) 



(9) 
(10) 

(11) 
(12) 

(13) 
(14) 



B. Global symmetries and conserved quantities 



The action |T]) is invariant under the following three-parametric group of global transformations 



ii)->toe 2 ( 5 + A ) ; p ^pc s , r ^ rc x + v, a^cre x , 



-26 



They generate three conserved Noether currents 

/ dL 



J > := {d¥^' A - L ) d ° r 



9=0 



, 9 r Jg = 0, 



5=0 



(15) 



(16) 



where 4> stands a,w,p,<p, /, and g = 8, u, A. The conserved quantity corresponding to the parameter v is the 
Hamiltonian 



1 



1 



H = ^[p'(wpY 1] + a'wpp' - 2aaa- V 2 w)'(3V 2 - l)4>'^ - ^ W p 2 <j>' 2 + I + ~q 2 m e 2a t (17) 

Z Z Z O Z p 



This is known to vanish on shell for diffeomorphism invariant theories, H = 0. The Noether current corresponding to 
the parameter 5 leads to the conservation equation J$ — 0, where 



awp 2 (j)' (a 2 w)'p 2 



J 8 = 



V -^- + 2qJ + 2a(ff- 1 {(wp' 2 - 1) [(a 2 w)' ~ 2aa 2 wc/)'] + 2aw(a 2 w)' pp' >'} e 2 <^) , (18) 



which is an Abelian counterpart of the equation given in [60j . The value of this integral depends on solutions. The 
third integral corresponding to A is trivial: J\ = —rH, its existence implies H = 0. 



1 We use this occasion to correct a misprint in Ref. [60ll : the factor e 2 * is missing at the right hand side of the Eq. (23). 
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The above integrals of motion allows one to reduce the order of the system by two. Fixing the gauge (e.g. a = 1) 
one has three second order equations for w, p, 4> with q e , q m entering as parameters of this six-order system. Using 
the integrals, the system order can reduced to four, with one parameter more (the fixed value of Js). For q e = one 
can further reduce the order to three. Introducing, for instance, new variables 

exp(u>), p -> exp - ^) , (j)^ <j)- —W, (19) 

we exclude from the system the variable w (while w' and w" still persist). For numerical computations we use the 
initial six-dimensional system, checking the constancy of the integrals of motion to control accuracy of the calculation. 

The space of solutions is invariant under a four-parameteric group of global transformations which consists in 
rescaling of the electric charge 

q e ^q e e 2S , q m ^q mi (20) 
(leaving the magnetic charge invariant), rescaling and shift of an independent variable 

r _> rc^ +s +v, (21) 
and the following transformation of the field functions: 

w^wef, p^pe 5 , a^ae x , ->•<£+-, f->fe$- s+x . (22) 

a 

Note that the Lagrangian is rescaled under this 4-parameter transformations as L — > e x L, so the action |T]) remains 
invariant provided 

p = -2(6 + A), (23) 

in which case we go back to (|15p . The shift v is trivial and the symmetry related to A can be frozen by the gauge 
choice cr=l. Therefore, physically interesting transformations are generated by p and S forming the subgroup which 
we denote as G(p, S). 

C. Turning points of p(r) and the gauge choice 

Re-paramctrization of the radial variable r allows to eliminate one of the three metric functions. There are two 
convenient gauge choices: the Schwarzschild gauge p — f, in which the radial variable is the radius of two-spheres: 



w(r 

and the GHS gauge [59| a = 1: 



ds 2 = -a 2 {f)w(f) dt 2 + -rj^ + f 2 dflj, (24) 



dr 2 

ds 2 = -w(r)dt 2 + —— + p 2 (r)dQ 2 . (25) 
w(r) w 2 y ' 



The coordinate transformation relating these gauges reads: 



1 dp(r) 

f = p(r), a 2 (f)w(f) =w(r), —-7 = —, — ■ (26) 

cr(r) dr 

It becomes singular at the turning points of the function p(r) where the derivative p'(r) — 0, so solutions containing 
such turning points can not be described globally in the Schwarzschild gauge. We will see later on that the presence 
of turning points is typical for the DGB system, so the GHS gauge is preferable. 

D. Dilaton black holes and the GB term 

In the theory without the GB term, a = 0, an electrically charged asymptotically flat black hole solution for an 
arbitrary dilaton coupling a reads (in the gauge (7 = 1) [59| : 

W[ r) = (l - T±) (l - r -) & , P (r) =r(l- r -)^, e^ = e 2 ^ (l - ^)"^ . (27) 
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The mass and the electric charge are given by 



""T + T^T- «*-«**" jt^. (28) 



For a = this reduces to the Reissner- Nordstrom solution, which in the extremal limit r + = r_ = th , 

, 2 , o dr 2 



fl-^ dt 2 + ^ 2 +r 2 rfO|, (29) 



has a degenerate event horizon AdS2 x ,S 2 . Note that for a = the GB term decouples from the system, so this 
solution remains true in the full theory with a^O. 
For a ^ 0, the extremal limit r + = r_ = r# is: 

d S 2 --(l-^) T ^^ 2 +(l-^-) _T ^ dr 2 +r 2 (l-^)^df! 2 . (30) 

At r = Tu the radius of the two-spheres shrinks, so we have a point-like singularity. The Ricci scalar in the vicinity 
of this point diverges near the horizon r = rjj- 

R ~ (l + a*)*T*V-T) ' (31) 
as well as the dilaton exponential for an electric solution: 



„2ad> 2. 

e v = e 



— 2a , 

a^o„ (l - ^ . (32) 



Substituting the general dilatonic black hole solution (|27p to the GB term we obtain the following value at the 
horizon r = r + 

e 2a *£ GB \ r=r+ ~ (r+ - r^)- 2 ^! 11 . (33) 

This expression diverges in the extremal limit r + — > r_. Thus, it is not possible to treat the GB term perturbatively 
expanding in a in the vicinity of the extreme dilaton black hole. In other words, one can expect that the GB term 
will substantially modify the dilaton black hole solution in the extremal limit. 

Summarizing the above information, we see that the DGB gravity admits the black hole solution with the degenerate 
event horizon of the AdS2 x S 2 type for a — (the Reissner- Nordstrom solution (|2"9"]) , and does not admit the extremal 
dilaton black hole (|30[) for a =/= even as an approximation for small a. So, the intriguing question arises, whether 
the branch of degenerate black holes exists in the DGB gravity which starts at a = in the parameter space and 
continues to non-zero a. In the next section we analyze this possibility in detail both analytically and numerically. 



III. DGB BLACK HOLES WITH AdS 2 x S 2 HORIZON 

We are looking for asymptotically flat solutions in the DGB theory for which the metric function w{r) has double 
zero at some point r = th and does not have singularities for r > rjj. To attack this problem numerically, one has 
to prove first that such solutions exist locally in the vicinity of the horizon r = th- We will show that this is true, 
provided some restriction on the parameters is satisfied. 



A. Near horizon expansion 



Assuming the GHS gauge a = 1, consider the series expansions around some point r = Th (supposed to be a 
horizon) in powers of x = r — r# : 

oo oo oo 

w{r)=Y,w k x k , p(r)=J2pkX k , P(r) := e 2a *^ = ]T P k x k . (34) 

fe=l k=Q k=Q 
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The function w starts with the linear term (vanishing of wq means that r — rjj is a horizon), two other functions 
have general Taylor's expansions. 

Substituting these expansions into the equations of motion ((9lfT2|) we find local solutions of two types. The first 
type solution has W\ ^ 0, i.e. the function w has simple zero at r — rjj- This corresponds to a non-degenerate horizon 
of the Schwarzschild type. Such local solutions and their numerical continuation to infinity were considered in some 
particular cases in refs. (U H, EH, |46|. Here we give more general expansion valid for both electric and magnetic 
charges present (in the gauge a = 1): 

T 2a 2 [a(qt + P 2 q 2 J + (4rP a 2 + p 4 )q 2 + P^P a 2 - p 4 + Aaq 2 )q 2 m ] + Tp® ^ ^ 

W{r) = tiPi 12p^P 2 a X +0{X } ' 

M , (Po-Po - - g| - 2aP T)P 1 2 

P{r) = P + Pix + O(x 2 ), (35) 

where pa, Po and Pi are free parameters and T satisfies the equation 

48a 3 a 2 P 2 T 2 + {pt - 16P aV [3p 2 P - 2(q 2 + P 2 q 2 m )} } V (36) 
+ 2a 2 [2a(qj + P \ 4 J + pUpI - l2P a a)q 2 e - p 2 P 2 {p 2 + l2P Q a)q 2 m + ^P 2 q 2 e q 2 n + 6P 2 ap 4 ] = . 

This quadratic equation has two roots T = T± depending on parameters q e ,q m , POi Pot P- The above local solutions 
exists for such values of parameters for which r ^ 0. 

The second class of local solutions has wi = 0. Vanishing of w\ means that the horizon is degenerate. Such an 
expansion contains only one free parameter with fixed charges. This family is disconnected from the family (I35| and 
it was not noticed so far: 



w ( r ) = X ~ - lAm [3(« 2 - !)<4 + 6«(3a 2 - 2)q 2 m + 4a 2 (5a 2 - 3)] x 3 + 0(x 4 ), 
Pq baa z pQ 



p{r) = Pa + [(a 2 - l)q 4 m + 2a(3a 2 - 2)q 2 m + W(a 2 - 1)] x + 0(x 2 ), 

P(r) = ^IJ+^ + ^ 2 )- (37) 

Here q m is the magnetic charge, Pi is a free parameter, and po is the physical radius of the horizon depending on 
charges as follows 

Pi = 2 «f« + f. (38) 

Note that the dilaton coupling constant enters these expansion only through a 2 , so the space of solutions is symmetric 
under a — > —a, <j) ~ * ~4>- I n what follows we will assume a > 0. 
The values of the integrals of motion corresponding to ([37j) are 

H = w 2 [q2 - Po + ^ - p2] ' Po = w^Kr m 

Js = 2g e / , (40) 

where /o is the value of the electrostatic potential on the horizon. 

From our previous analysis of the Einstein-Maxwell-Dilaton (EMD) black holes and the above form of the local 
solution one can draw the following conclusions: 

1. Black holes with A1IS2 x S 2 event horizon do not exist in the a ^ EMD theory without curvature corrections 
(a = 0). 

2. Such solutions do not exist in absence of the electric charge, while the presence of the magnetic charge is optional. 
S-duality is thus broken as expected. 

3. This local solution is not generic (the number of free parameters is less that the degree of the system of differential 
equations) . 
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For simplicity, in this paper we will focus on purely electric black holes. In this case the expansions simplify and 
we can give some further terms: 



w(r) 
p(r) 
aP{r) 



1 

7 



Po 



2(5a 2 -3) / aPi 



a 2 pl 



173a 6 - 269a 4 + 99a 2 - 27 / aPi 



3(5a 2 - 3) 



+ 0(x 5 



1 + (a 2 - 1) 



aPi 
o?p 2 



2a 2 (a 4 -6) ( aP l 
(5a 2 - 3) 



a 2 p 2 



= Pi 



4 +a 



aPi 

~ 2 2 

a Po 



a 2 (a 4 -5a 2 -3) ( aP l 



(5a 2 - 3) 



a?pl 



+ 0(x 3 ), 
+ 0{x 3 ). 



The electric charge is related to the horizon radius as 



Qe 



2V^' 



(41) 



(42) 



We can arrange higher order terms in such a way, that a 2 enters in denominator only in powers of the combination 
Pi / a 2 . This facilitates taking the limit a — > 0. We know that in this limit there exists an exact solution which has 
the near-horizon expansion of the type (141|) , namely the extremal Reissner-Nordstrom solution. Thus we expect that 
for the asymptotically flat solutions 



Pi -» 0, 



0. 



(43) 



Numerical calculations show that this is indeed the case (see Sec. C). 

Another subtlety is related to the limit of GB decoupling. Obviously, our expansion fails in this limit: for a finite 
charge q e one has p\ — 2q e y/a — > and, consequently, the expansion coefficients will diverge. The reason is that our 
expansion for w is incompatible with that for the dilaton black hole of the Einstein-Maxwell-Dilaton theory. This 
reflects again the absence of the black hole solutions with the AdS2 X S 2 horizon in the theory without curvature 
corrections. Substituting q e defined by ([42]) into the equations of motion (j9l[l~2"|) . one can see that the GB coupling- 
parameter a enters always in the combination ae 2a ^ r \ Thus, shifting the dilaton is equivalent to rescaling the GB 
term (this was in fact clear already from the action |T])). Note that in the case a 2 = 1 the linear term in the expansion 
of p vanishes, pi = 0, implying that there is no regular transformation to the gauge p — f. 

Therefore, a purely electric local solution with a fixed value of charge q e contains one free parameter: the dilaton 
derivative Pi. An important issue is to determine the correct sign of Pi. To be able to interpret the region r > Th 
as an exterior of the black hole, one has to ensure positiveness of the derivative p' at the horizon. From the above 
expansion one finds 



'I (a 2 - l)aPi 
p \x=o = 5 > 0. 



(44) 



Thus, we should take positive Pi for a 2 > 1 and negative Pi for a 2 < 1. It is convenient to introduce the sign 
parameter which ensures this: 



C = 



_P_ 

\Pi\ 



a 2 -l 



Now define the following combination of po and Pi : 



b = 



<*\Pi\ 
a 2 p 2 ' 



(45) 



(46) 



It is easy to see, that free parameters enter the series expansions near the horizon only through this quantity. Consider 
now the transformations of the expansion parameters under symmetries of the solution space (|21I22[) . Since Pi is the 
first derivative of the dilaton exponent, one finds that under (^-transformation 



Pi 



Pie fl 



be 



(47) 



so the quantity bx remains invariant. Thus the full set of local solutions can be generated from one particular 
solution with po — 1, Pi = 1, which we will call the normalized local solution, by the symmetry transformations with 
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5 = — lnpo and p = 2 \n(bpo) from (|22| . i.e. by the group element G(2 ln(bpo), — In po). The normalized local solution 
does not contain free parameters at all: 

2 2(5a 2 -3) 3 (173a 6 - 269a 4 + 99a 2 - 27) 4 ^. 4 . 
p(r) = 1 + ,(a 2 - l)x - 2 ° 5 X~3) ) x2 + < 49 > 



1 . , (a 4 -5a 2 -3) , 

4 + ^-^ )+ (5a 2 -3) 



aP(r) = 1+q (r-r H )+ v /c 2 V + Ofr 3 ). (50) 



Note the presence of the sign function «; in the odd power terms. The electric charge corresponding to the normalized 
local solution is q e = l/(2y/a). 

B. Asymptotic flatness 

We are looking for asymptotically flat global solutions which satisfy the conditions w — ► 1, p/r — ► 1, <\> — > const 
as r — > oo. The subleading terms should be expandable in the power series of 1/r. The local solution with these 
properties turns out to be three-parametric, depending on the ADM mass M , the dilaton charge D, and the asymptotic 
value of the dilaton : 

. 2Af aQ 2 _ 3 

u>(r) = 1 1 p- + 0(r J ), 

r 

L> 2 D(2MD - aaQl) 
p(r) = ^ - - — ^ — + °( r )• ( 51 ) 

, , . , D 2DM-aaQl _ 3 , 

Wh6re Q. = ff»e-*-. (52) 

The dilaton charge can be also read from the asymptotic expansion of the dilaton exponential: 

= 1 + 2 -f + ™>W + M)-™ 2 <X + 0(r - 3) . (53) 

The asymptotic values of two integrals of motion are: 

H = \ ( Woo p'l - 1) , (54) 

Js = 2 q e foe — M ——. (55) 
a 

Behavior of the global solution which starts with the normalized local solution (|48I49I50[) at the horizon depends 
only on the dilaton coupling constant a. Its existence for all a is not guaranteed a priori. But, for some sufficiently 
small values of a, we find numerically that all three functions vary smoothly with increasing x, so that w and the 
derivative p' stabilize at infinity on some constant values u>oo ^ 1, p' x ^ 1. Then, using the symmetries (|21I22[) 
of the solution space, one can rescale the whole solution to achieve the desired unit values for these parameters. 
More precisely, the relevant subgroup of rescalings is two-parametric: G(p,5). As we have argued, two parameters 
p,S effectively replace the parameters po,Pi of the (non-normalized) local solution (|4"Tj) . So one could expect that 
rescaling of the solution so that w^o = 1, = 1 would fix both quantities po,Pi on the horizon. But from the 
Hamiltonian equation H = with H given by the Eq. ()54|1 it is easy to see that one must have Waop 1 ^ = 1 for any 
solution such that w — > u>oo , p' — > p'^ asymptotically. Therefore it is enough to perform one but not two independent 
rescalings in order to get Woo = 1, p'^ = 1. Indeed, under G{p, 5) 

w^wc^, po^poe 5 , Pi^Pie^/ 2 , wp' 2 ^wp' 2 . (56) 

Since the choice of p, S is equivalent to the choice of po, Pi, an invariance of the product wp' 2 under G{p, 8) means that 
the solution starting on the horizon with any po,Pi will reach at infinity the values Woo, p'oo satisfying Woop'^ = 1. 
Therefore, taking p = — In Woo, we will achieve simultaneously w m = 1 and p M = 1. This means that asymptotically 
flat solutions still form a one-parameter family, a parameter being the electric charge q e . 
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C. Numerical analysis 

Since we know that the desired global solution exists for a — 0, we start with the local solution at the horizon 
with small a and look for numerical solutions which fit the asymptotic expansions (|51[) . For sufficiently small a global 
solutions exist indeed, and, as we have explained, two basic conditions at infinity w — 1, p' = 1 fix only one of the two 
parameters po and P\ at the horizon. It will be convenient to leave po (defining the electric charge) arbitrary, and to 
fix Pi. We will also choose the value of the GB coupling a = 1. Then the black hole mass can be found numerically 
from the asymptotic expansions (|5ip together with the dilaton charge and the asymptotic value of the dilaton. 

Typical coordinate dependence of the metric functions and the dilaton exponent are shown in Fig. Q]for some values 
of the dilaton coupling a. Solutions exist for 

< a < a CI , a cr ~ 0.488219703. (57) 

Let us discuss in more detail behavior of solutions at the ends of this interval. As expected, the parameter b of the 




° 20 40 60 80 100 120 140 0£ 1 1^5 2 1 2 3 4 5 6 



FIG. 1: The functions w(x), p'(x), P(x), x = r — th, for po = 1 and some values of the dilaton coupling constant: a — 0.1 
(Pi = -0.01) — thin line, a = 0.4 (Pi = -0.446) — normal line, a = 0.45 (Pi = -1.18) — thick line. 

local solution at the horizon tends to unity when a — > 0, as shown in Fig. 12.11 This means that the first Taylor 
coefficient in the expansion of p{x) becomes equal to unity (note that the sign function c = —1 as a — > 0), while all 
higher coefficients are zero. Therefore, assuming po = th, we find that p = r globally. Similarly, all terms in the 
expansion of P{x) vanish in the limit a — > except the constant Po, so the dilaton exponential tends to the constant 
value P — Pq/4. Correspondingly, we find 

lime 2a ^/q e -> 1/2. (58) 

For w all Taylor's coefficients in (|4ip are non-zero and the whole series exactly reproduces an expansion 

w(r) = (l-^y =z 2 -2z 3 + 3z 4 + 0(z 5 ), z=(r-p )/p . (59) 

Thus, our family of solutions begins with the extremal Reissner-Nordstrom metric for zero dilaton coupling a. 

With fixed horizon radius po, the mass and the dilaton charge of the black holes increase with the growing dilaton 
coupling constant tending to infinity when a approaches a CI . The dilaton exponent, on the contrary, tends to zero in 
this limit. Using the symmetry of the solution space under (5-transformation, one can generate the sequence of solutions 
with different electric charges q e and correspondingly with different masses, dilaton charges and ^oo. Since variation 
of the electric charge is essentially equivalent to variation of the unique parameter po in the horizon expansion, it is 
clear, that using ^-transformation we will generate all extremal solutions. Under this transformation the mass and 
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FIG. 2.1: Dependence of the parameter b on the dilaton coupling constant a: b tends to unity for a — > ensuring continuous 
transition to the extremal Reissner-Nordstrom solution. 

FIG. 2.2: Numerical curves k,M{a) = M 2 /q e , fczj(a) = D 2 /q e , k^(a) — e 2a '* 30 /g e for po = 1- The mass curve starts with the 
Reissner-Nordstrom value for a = and diverges as a — + a cr . The dilaton charge increases from zero to infinity, while the 
dilaton exponential monotonically varies from the value 1/2 at a — to zero for a — » a CI . 



FIG. 2.3: The quantities — ry- and — — -p. as functions of a: both tend to unity at the ends of the allowed interval of a. 



the dilaton charge scale as e 5 , while the electric charge and the dilaton exponent e 2a(?io0 scale as e 2S . Therefore the 
ratios 

M 2 D 2 e 2a <t>°° , . 

km = , k D = — , ks = (60) 

q e q e q e 



depend only on a. Their numerical graphs are presented on Fig l2.2l 

As we already discussed, the metric for a = is known analytically. For a in the vicinity of a cr the analytic solution 
is not known, but one finds that the behavior of the mass, the dilaton charge and the rescaled electric charge Q e 
exhibits some similarities at both ends of the allowed interval of a. Namely, the following two ratios stabilize at unity 
for both a — > and a — > a cr (Fig. 



aM 



D 

This corresponds to fulfillment of the following condition 



VI + a 2 M 



1 (61) 



a 2 M 2 + D 2 = -^Q 2 e , (62) 
1 + a z 

which is reminiscent of the BPS condition for electrically charged black holes of the EMD theory. This feature is 
similar to that in another stringy generalization of EMD theory in which the Maxwell action is replaced by the 
Born-Infeld action, but no GB term is introduced (6l| . 

For the values of a outside the allowed interval, solutions starting with the AdSi x S 2 horizon are not asymptotically 
flat, but singular. For them the metric function p(r) has a turning point at some finite radial coordinate r = r t , such 
that p'(r t ) = 0, p"{r t ) < (Fig. [3]). This point is regular, but at a finite proper distance from it one encounters the 
singular turning point r s , such that p'{r s ) > 0, p"{r 8 ) = oo, where all variables have a square-root singularity, being 
expandable in terms of V r ~ r s- 

w = w s +w 1 y + w 3 / 2 y 3/2 + 0(y 2 ), y = r-r s , 

P = Ps + Piy + P3/2y a/2 + 0(y 2 ), (63) 
P = P s +P iy + P 3/2 y 3 / 2 + 0(y 2 ). 
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9 0.2 0.4 0.6 0.8 1 1.1 1.2 1.3 1.4 1.5 



FIG. 3: The metric functions p(x),w(p), (x = r — ru) for the value of the dilaton coupling constant a > a cr . The function p 
has a turning point r — rt after which the naked singularity is met at a finite affme distance (r = r s ). Dotted lines correspond 
to another solution branch which can be matched in the singularity. Numerical curves are presented for a = 0.5 and po = 1, 
the corresponding value of Pi being Pi = —7.746. 



This local expansion contains four free parameters w s , p Sl pi, P s , while other coefficients read: 
Aa 2 P s ( 1P 2 s P s +q 2 e )- A 2 w s pj (Ps - 4A Pl w s )p 3/2 

Wl = AU2 2 2 r (R\ \ TaT ' Pl = A ' P 3/2 = o » ( 64 ) 

AP£a z pj [pi(QAw s pi - Ps) - 2Aj 2j 
{2Apj [2A Pl w s ( 7 + 4) - Ps (5 7 + 4)] - 16a 2 P sPll ( 1Pl P s p 2 s + g 2 e ) + p e sPl } w B p 3/2 
Wm 8p^p sV(2A(37 + 2) _ psPl) > ^) 

where 7 = p\w s — 1 and A satisfies the equation: 

A 3 {8w sP 4 s [«w>?(157 + 9) - 1] } - A 2 24 Ws p5 pi(37 + 2) + psf>i [ q ^ 2a 2 P sl + P 2 s (48a 2 P 2 7 2 ~ p 4 s )] (66) 
+ A {2p 2 [p 4 s ( 7 + 6) + QGplws^P 2 ^] - ge 2 32a 2 P s (3 7 + 4) 7 } = 0. 

An expression for ,03/2 is too big and is not given here. Since the second derivatives are divergent at y = 0, the 
Riemann tensor diverges as well. The divergency is localized on a sphere of finite radius, and it is rather mild: Ricci 
and Kretchmann scalars behave as 

„ 3(Ap 3/2 w 3 + w 3/2Ps ) 1 uQ0jS 9(8p 2 /2 w 2 + p 2 s w 2 3/2 ) 1 

H -. 1=, H a/ 3jsn H ' ~ — -2 . (67) 

4p s y/V iGpj v 

The radial coordinate stops at r = r s , but using an appropriate desingularization of the system (see Appendix), one 
can glue another patch of radial coordinate r' S (r s , r/) at this point, extending the manifold through the singularity. 
This extension is shown by dotted lines in Figs. [3l It terminates at the final singularity r/. This situation is very 
similar to that described in ref. [46j for an interior region of the non-extremal DGB black hole. 



D. Thermodynamics 



The temperature of the extremal DGB black hole is zero, as for the extremal solution without the GB term: 



T = — 
2tt 



dr 



1 



2tt P 2 



(r - r H )\ r=rH = 0. 



(68) 
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To calculate the entropy we apply Sen's formula [33, Ho| appealing to the near horizon data. Using (|4Tj) we can write 
the near horizon solution as 



PO ( r - r H) 2 



ds * = ^di A + - ru dr 2 + pgdfl 2 ., 



fln|°, F [2] = ?#diAdr. (69) 
2a 4a Pq 



To apply Sen's formula [39j, [40(, we rewrite it as follows 

ds 2 = vi -r 2 di 2 + — J + U2^|, 
^» = u, F H = e, (70) 

where f — r — rH, t = t/p 2 . 

Now introduce the surface integrated Lagrangian density 

f(u,vi,v 2 , e) = / d8d<f)yfg L, (71) 

and evaluate it using the near horizon data (170|) 

/(«, «i,«a, e) = ~ (vi - va - 4ae 2au + e 2 ^e 2a ^ . (72) 

The entropy function _F is the Legendre transform of this function with respect to e: 

df =e V2 e 2au F = 2ir[qe-f(u,v 1 ,V2,e)}=ir(v 2 -v 1 +4ae 2au + e 2 —e 2a A, (73) 



or, in terms of g: 



F = Tt(v 2 -v l +Aae 2au + q 2 ^-e- 2au \ . (74) 

The entropy of the extremal black hole is given by the value of the entropy function F at extremality: 

dF dF dF 

ou av\ ov 2 



In our case, the extremality conditions (|75p read 

„ — 2au 2 ri 2 , „— 2au 2 r* ^ — 2au 2 , a /2au n 

e g — «2 = 0, — u 2 + e g Ui = 0, — e q v 1 + 4ae W2 = 0, (7o) 
leading to the solution 

v 1= 2^q, v 2 = 2^q, e 2au = -^=, (g>0). (77) 

2v« 

Comparing with the local solution in our previous notation, we get q = q e . Finally, substituting (|77|) in F one obtains 
the entropy 

S = 4nry/a q e = 2%(%, (78) 
which is precisely twice the Bekenstein- Hawking value. This is similar to the result of refs. [20l. [2ll. |24T] . 

IV. DISCUSSION 



In this paper, we have shown that in addition to charged black holes with non-degenerate horizons, the DBG 
four-dimensional gravity admits black hole solutions with the horizons of the AdS2 X S 2 type. These solutions form a 
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one-parameter family and exist in a finite range of the dilaton coupling constant a. New family of solutions branch is 
disconnected from the branch of non-extremal black holes which was studied earlier. Rather, it pinches off from the 
extremal Reissner- Nordstrom black hole which is a solution of the full EMDGB theory for a = 0. Starting with zero a, 
we were able to find global black hole solutions interpolating between AdS2 x S 2 at the horizon and Minkowski vacuum 
at infinity for a below some critical value which was found numerically up to several decimals as a cr ~ 0.488219703. 
Near the critical value a — > a cr , the mass and the dilaton charge grow up, while their ratio saturates the BPS bound 
of the EMD black holes. Similar feature was observed for the charged black holes in the Einstein-Born-Infeld-dilaton 
(EBID) theory [H[. 

It is worth noting that the family of electrically charged extremal black holes in the EMDGB theory is one- 
parametric (<7 e ), while the family of the corresponding extremal solutions in the EMD theory is two-parametric (with 
the parameters q e and <^oo). An asymptotic value of the dilaton is no more a free parameter when the Gauss-Bonnet 
term is included, moreover, the dilaton exponent e 2a ^°° at the threshold a = a CI tends to zero for any finite value of 
the charge q e . Therefore, modification of the extremal dilaton black hole by higher curvature term consists not only 
in stretching its horizon to a finite radius, but also in fixing the value of the dilaton at infinity. 

Our model can be viewed as a truncated heterotic string effective theory in four dimensions. Whilst it does not 
include all quadratic curvature terms, it still exhibits features relevant to more complete models, in particular, it 
predicts correct entropy for extremal black holes which is twice the Bekenstein-Hawking entropy. The dilaton varies 
from a finite value at the horizon to some different finite value at infinity. The existence of the threshold value of the 
dilaton coupling constant under which the global solutions cease to exist is an interesting new phenomenon which may 
be related to the string-black hole transition as described in [62] 2 . We think that our model as well as the EBID model 
[6l| (both incorporating typical stringy features) can be regarded as simple toy models describing the string-black 
hole transition. 
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Here we clarify the numerical procedure which allows us to continue solutions through the singular points. Rewrite 
the system (0 - [T2")) as a matrix equation of the first order 



where X is the six-dimensional vector consisting of the primary dynamical variables w, p, e a * and their first derivatives 
with respect to the radial coordinate. The system (IA1|1 has a regular solution provided det A ^ 0. When the solution 
approaches some point r s where det A — > 0, the derivative X' diverges as 0(l/det A). In order to continue the solution 
through this point, we choose a new independent variable a satisfying the condition 
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APPENDIX A: DESINGULARIZATION AT THE TURNING POINT 




(Al) 



da 

Then in terms of a the matrix equation (|A1|) can be rewritten in the regular form 



r(a) — — — cx det A. 



(A2) 



AX - Br = 0. 



(A3) 



2 We thank Miguel Costa for bringing the paper [62$ to our attention. 
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This desingularization is achieved by extending the set of unknown functions to seven, considering the radial coordinate 
as a function r{a). Denoting the seven-vector (X(a), r(a)) as Y(<r), one can see that the tangent vector has the unit 
Euclidean metric norm: 



dY 



da 



1, 



(A4) 



provided the equation (|A3j) holds 



|detA| (det.4 2 + (X'detA) 2 



(A5) 



Using this desingularization one can continue the solution through the singular turning point. This procedure is 
similar to one used in [631 ] . Geometrically this means gluing another coordinate patch to the solution at singularity. 
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